Rydberg excitation of Bose-Einstein condensates 



Rolf Heidemann0 Ulrich Raitzsch, Vera Bendkowsky, Bjorn Butscher, Robert Low, and Tilman PfaiJll 
5. Physikalisches Institut, Universitdt Stuttgart, Pfaffenwaldnng 57, 70569 Stuttgart, Germany 

(Dated: 29 October 2007) 

Rydberg atoms provide a wide range of possibilities to tailor interactions in a quantum gas. Here 
we report on Rydberg excitation of Bose-Einstein condensed *^Rb atoms. The Rydberg fraction 
was investigated for various excitation times and temperatures above and below the condensation 
temperature. The excitation is locally blocked by the van der Waals interaction between Rydberg 
atoms to a density-dependent limit. Therefore the abrupt change of the thermal atomic density 
distribution to the characteristic bimodal distribution upon condensation could be observed in the 
Rydberg fraction. The observed features are reproduced by a simulation based on local collective 
Rydberg excitations. 



The ability to tailor interactions among atoms in a 
quantum gas is one of the main strengths of ultracold 
atoms as model systems for other research fields such as 
condensed matter physics. So far mostly effective contact 
interactions have been used which can be changed over a 
wide range in strength and even in sign. Two years ago, 
the long-range magnetic dipole-dipole interaction was ob- 
served between atoms in a ^^Cr-BEC [l|. It waspro- 
posed to induce dynamic dipole-dipole interactions [4I or 
a 1/r interaction potential Q between ultracold atoms 
by light. To induce a static electric dipole-dipole interac- 
tion, it was suggested to admix a Rydberg state with a 
static electric dipole moment to the atomic ground state 
by photo-excitation Q. 

Rydberg atoms and ground state atoms are predicted 
to form unusual weakly bound molecular states [1, 0]. 
As the average interparticle distances in BECs are of the 
same order of magnitude as the molecular binding dis- 
tances, a BEC seems to be the ideal starting point for 
photoassociation of these molecules. 

A Rydberg atom also constitutes an impurity in the 
BEC related to Ref. which could be manipulated by 
electric fields. The interaction could lead to an agglomer- 
ation of ground state atoms around the impurity as was 
calculated for ionic impurities Q. 

As the nature and strength of the interaction among 
Rydberg atoms can be tailored by electric fields and mi- 
crowave fields e.g. from van der Waals type to dipole 
dipole type or an isotropic potential they 
provide a new tool for many body quantum physics. 

In this Letter we describe the excitation of atoms to a 
Rydberg state while undergoing a phase transition from 
a thermal gas to a Bose-Eiirstein condensate. We present 
a model based on collective excited states to simulate 
the observed Rydberg fraction across the phase transition 
which agrees qualitatively with the observed data. 

Using the setup described in , we magnetically trap 
^"^Rb-atoms in the 5Si/2,^ = 2,mp — 2 state and pro- 
duce samples from thermal clouds to Bose-Einstein coir- 
densates by means of forced RF-evaporation. We vary 
the temperature from 5 /uK down to 200 nK crossing the 
condensation temperature at around 700 nK. After this 



preparation, the atoms are subject to a two-photon Ryd- 
berg excitation via the 5P3/2 state to the 43Si/2 state. 
The duration of the square pulses of excitation light was 
varied in this experiment between 170 ns and 2 /is. To 
avoid significant absorption and heating due to sponta- 
neous photon scattering, the light is blue detuned by A= 
483 MHz from the SFg/a,^ = 3 level. Thus only one 
photon per 100 atoms is scattered for the longest ex- 
citation time. For the following experiments the Rabi 
frequency fii on the 5S-5P transition is 11 MHz. A Rabi 
frequency ^2 for the 5P-43S transition of 9.7 MHz ^ 
results in a two-photon Rabi frequency f7ii72/(2A) of 
110 kHz. The thermal motion of the atoms is negligible 
on the /^s-time scale of the experiments, but interactions 
among the Rydberg atoms can lead to collisions and ion- 
isation. By choosing a Rydberg state with repulsive in- 
teractions and applying an electric extraction field, the 
unwanted effects of ions are avoided [l^ . 

The Rydberg excitation is followed immediately by a 
field-ionisation of the Rydberg atoms and their detection 
on a microchannel plate. Up to three such sequences of 
excitation and detection are applied to one sample fol- 
lowed by absorption imaging of the remaining ground 
state atoms. From these calibrated detection methods, 
the total atom number Ng and the fraction / = Nfi/Ng 
of atoms in the Rydberg state are derived. This normal- 
isation of the Rydberg atom number Nf; is useful since 
the total atom number drops significantly from 10^ to 
10^ while reducing the temperature. The temperature of 
each sample is derived from the size on the absorption im- 
age after an expansion during 20 ms time of fiight. The 
thermal clouds are fitted using a Bose-distribution ne- 
glecting interaction among the ground state atoms. The 
BEC component is fitted with a Thomas-Fermi distri- 
bution. From the parameters of the fit and the known 
trapping potential, the in-trap density distributions are 
calculated and used for simulating the Rydberg excita- 
tion. The peak density and the condensate fraction are 
plotted versus temperature in Fig.[T] where no effect of 
the different excitation times oir the atomic distribution 
is visible. The bimodality of the density distribution will 
prove to be crucial for the understanding of the results 



2 



(a) 



(a) 



(b) 



full distribution 
thermal fraction 





10'° 


(b) 






0.6 


tion 


0.4 


u 
nl 




te fr 


0.3 


;nsa 


0.2 


■D 
c 




o 
o 


0.1 



T = 1x 1970 ns 
T = 2x 670 ns 
T = 3x 370 ns 
T = 3x 320 ns 
T = 1x 170 ns 





0.2 



).5 1 2 5 

temperature (^iK) 



FIG. 1: (a) Dependence of the peak density of ground state 
atoms on the temperature. The crosses ( x ) indicate the peak 
density of the thermal component, the dots (■) the peak den- 
sity of the bimodal distribution (where a BEG component 
exists). In (b) the condensate fraction is plotted versus the 
temperature. The colors indicate the length of the Rydberg 
excitation pulse, which was applied in trap before taking the 
absorption picture after a time of flight of 20 ms. The num- 
ber of identical excitations within one sample is denoted by 
the multiplier in the legend. Note that no systematic effect of 
different excitation times on the distributions is visible. The 
same samples delivered the data presented in Fig. [2ja). 



of the Rydberg excitation. When the temperature is de- 
creased below the critical temperature, the central den- 
sity increases abruptly by a factor of about 4. When 
the temperature is reduced further, more atoms are con- 
densed into the small constricted region of the BEC re- 
ducing the density of the thermal component. In the con- 
densate the average atomic distance to the next neigh- 
bour is 150 nm 0.55ng ^l?]) which is calculated from 
the peak atomic density Ug. 

In a recent work, we investigated the density depen- 
dence of the Rydberg excitation dynamics in magnet- 
ically trapped clouds above the condensation tempera- 
ture The Rydberg atom number as a function of 
excitation time initially increases linearly with a slope R 
and saturates to a value A'sat after a time Ts = Nsat/R- 
From this previous experimental observation we know 
that the saturation is reached faster with increasing den- 
sity and that the saturation Rydberg density depends 
only very weakly on the density of ground state atoms 
Tig. In this Letter we investigate the excitation dynam- 
ics across the condensation temperature. Due to the bi- 
modal density distribution in a partially condensed cloud 
the time scales and saturation Rydberg fractions are dif- 
ferent in the condensate and thermal component. Figure 
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FIG. 2: Rydberg fraction as a function of temperature for 
different excitation times. The data sets for each excitation 
time are color-coded and for clarity shifted such that the dot- 
ted horizontal lines each indicate 10~^. The measurements 
in (a) are smoothed for each excitation time giving the solid 
lines as a guide to the eye . The solid lines in (b) show the 
calculated Rydberg fractions based on the superatom model. 
The dashed lines show the result of the simulation taking only 
the thermal part of the atomic cloud. This confirms the Bose- 
Einstein condensation as the origin of the kink in the Rydberg 
fraction. 



O^a) shows the Rydberg fraction after different excita- 
tion times as a function of temperature. The lines are 
the result of a smoothing intended to guide the eye to 
the main features of the measurement: At medium ex- 
citation times (320 and 370 ns) the Rydberg fraction de- 
creases with proceeding formation of the condensate i.e. 
reduced temperature. For the shortest and longest ex- 
citation times, no such kink in the Rydberg fraction is 
visible at the critical temperature. 

The excitation times in this experiment are chosen such 
that for the shortest time (170 ns), both components are 
just slightly blocked but not saturated. For longer ex- 
citation times, the excitation is already saturated in the 
condensate component while it is not yet in the thermal 
part. For the longest excitation time (1970 ns), the ex- 
citation is saturated in both components. As a start the 
following intuitive picture may be used to explain the 
reduction of the total Rydberg fraction at medium exci- 
tation times: With increasing BEC fraction, more atoms 
are condensed into the BEC where the Rydberg fraction 
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is significantly lower than in the thermal cloud due to its 
higher density, this also lowers the total Rydberg frac- 
tion. This effect dominates over the smaller increase in 
the Rydberg fraction in the thermal cloud as its density 
decreases (see Fig.[TJa)). A deeper comprehension is ob- 
tained by the following model based on collective states. 

In this experiment the van der Waals interaction 
among the Rydberg atoms blocks further excitation 
within a blockade radius rj, around one Rydberg atom 
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17l |. This limitation to one excitation per block- 
ade sphere leads to a collective Rabi oscillation which 
is enhanced in frequency by a factor -x/iV. In this sense 
the atoms per blockade sphere act like a 'superatom' 
[lil with a -s/^V larger transition matrix- element. In 
the regime where the size of the sample is larger than 
the blockade radius rf,, quantum correlations intercon- 
nect the whole sample. Thus the system with Ng « 10^ 
atoms and Nn « 10'^ excitations is inaccessible to di- 
rect ab initio simulations. The model described below 
reduces the quantum correlations to spatial correlations 
and locally collective states by assuming Nfi independent 
superatoms which oscillate at their respective collective 
Rabi frequencies a/TV Qq ■ We will assume the superatoms 
to arrange in the close packing of a face-centered cubic 
structure with 12 next neighbours. In this hypotheti- 
cal lattice, the saturation density is nn = ^/2r^^ and 
the interaction energy at the center of one superatom is 
Z = 14.5 times the pairwise van der Waals ernergy. Z 
is slightly larger than 12 due to the atoms beyond the 
shell of next neighbours. In a local density approxima- 
tion the local number of ground state atoms forming one 
superatom is given by N{r) = ng{r)/nii{r) and density 
variations on the scale of rj, are neglected. The saturation 
density of Rydberg atoms can be derived from the block- 
ade condition which equates the van der Waals interac- 
tion energy with the collective Rabi frequency, where n 
is a constant expected to be on the order of one: 
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^ZCe-nlir)^Kh^/N{?jno. 



(1) 



The Ce-coefficient for the 43S'-state is -1.7 x 10^^ a.u. 
[l9| . As shown later, the blockade radius is in the mi- 
crometer regime, substantially larger than the mean in- 
teratomic distance. At these interatomic separations, the 
van der Waals interaction that is used in Eq. ([T]) is jus- 
tified. From this relation, the local saturation density 
nf{{r) and then the collective Rabi frequency Qc can be 
directly derived. 
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To illustrate the calculated bimodality in the satura- 
tion density of Rydberg atoms nii{r), Fig. [3] shows the 
density distributions as well as the distribution of the 
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FIG. 3: In (a) the radial cross sections of the density distribu- 
tions are shown for the ground state atoms 7ig(r) in red/light 
grey and for the simulated saturation density of Rydberg 
atoms n_R(r) in blue/dark grey. The cross section of the local 
Rydberg fraction /(r) is shown in (b). Each figure shows the 
calculated distributions for the lowest temperature of Fig. [5] 
and the dotted lines represent the thermal component without 
BEG. As can be seen, /(r) decreases with ng{r). Therefore in 
the BEG / is substantially lower than in the thermal cloud. 



fraction /(r) which shows a significant reduction at the 
position of the BEC component. The time dependence 
of the total number Nji{t) of Rydberg atoms can be ob- 
tained by integrating the oscillations of the superatoms 
over the sample: 
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The simulated behaviour of Nr(t) nicely reproduces the 
saturation curves shown in Ref. [13| . The derived scal- 
ing of Asat oc nn cx nj^ilfl^ with Rabi frequency is in 
excellent agreement with our previous observation where 
we found A^sat oc n° 07±0 02f]0-38±0 04[J|^ ^j^^ ^^^j^^g ^^^^ 

the peak density Ug is close. The calculated scalings of 
the initial slope of the excitation curve, R oc nn flc oc 
n^J^ffJ^, are in good agreement with the previous ob- 



servation of R 



,0.49±0.06^1.1±0.1 







[IJ. This simphfied 



model, despite its approximate character, explained well 
the observed scalings. Here we apply it for the situation 
of a changing density distribution ng{r) upon condensa- 
tion. The lines in Fig.l^Jb) show the calculated values 
for the Rydberg fraction using this model. The atomic 
density distributions were calculated from the informa- 
tion obtained from the absorption images and averaged 
over several pictures taken at different excitation times. 
The simulation uses two free parameters which is the 
propotionality factor k in Eq. ^ and the Rabi-frequency 
Oo- The factor k was adjusted to 0.3 such that the mea- 
surement is best reproduced for high temperatures. The 
Rabi frequency had to be reduced by a factor of 5.5 to 
reproduce the characteristic kink at the right excitation 
times. This reduction is reasonable as the frequency of 
the excitation lasers is only stable within 1.5 MHz be- 
tween adjacent excitations thus effectively reducing the 
Rabi-frequency in this model. Note that the simulation 
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was used beyond the large radius r, where the model 
gets unphysical due to nuir) > ng{r). For our parame- 
ters the effect on the total fraction / is negligible as only 
a fraction of at most 10~^ of the atoms are at these large 
distances where the collective behaviour ends and their 
single-atom excitation fraction is negligible. 

The scaling of the Rydberg fraction / = Nj^/Ng with 
density can be derived from the model for short, mod- 
erate and long excitation times: For short excitation 
times, when the dynamics cannot yet be distinguished 
from single-atom behaviour, the Rydberg atom number 
is proportional to Ng, therefore / is independent of Ug. 
For intermediate excitation times r, during the linear in- 
crease of Nji{t), f can be approximated by Rr/Ng oc 
tirVIcT /ug (X Ug . For long excitation times, in the 
saturation regime, the dependence of the fraction on the 
density is stronger: / cx N^^t/Ng (x nn/ug cx n'g^f(\ 
These considerations hold for a change in density with- 
out changing the shape of the distribution i.e. without 
crossing the critical temperature or for the thermal com- 
ponent alone. 

The simulation is able to reproduce the main features 
of the experimental data: The dependence of / on the 
density is increasing for increasing excitation time. For 
intermediate excitation times, / bends down below the 
condensation temperature. At long excitation times, this 
kink is covered by the increased dependence of / on Ug 
in the saturated thermal cloud. To emphasize the effect 
of the bimodal density distribution, the dashed lines in 
Fig.m^b) show the calculation which includes only the 
thermal component of the atomic cloud which does not 
reproduce the characteristic kink. According to this sim- 
ulation, the maximum density of Rydberg atoms in the 
BEC is 9xl0^^m~^, corresponding to a blockade radius 
of 5.4 /xm and 5 Rydberg excitations within the volume 
of the BEC. 

The experimental fraction / shows an overall increase 
with decreasing temperature which is not sufficiently re- 
produced by the simulation. We attribute this remaining 
discrepancy to a failure of the local density approxima- 
tion (LDA). The radial width of the gaussian distribution 
gets as narrow as 2 /im while the dynamics depends on 
a density that is averaged over the size of the blockade 
radius. For low temperatures this averaging effectively 
lowers the ground state densities and leads to higher Ryd- 
berg fractions than the simulation suggests. 

Note that we do not expect the change in the density- 
density correlation upon condensation to have a signif- 
icant effect on the described measurements as the de 
Broglie wavelength which is the length scale on which 
the bunching in the thermal cloud occurs is on the or- 
der of (< 350 nm) and therefore much smaller than the 
blockade radius. 

To conclude, we have demonstrated the Rydberg ex- 
citation of Bose-Einstein condensates and show that the 



momentum distribution of the condensate is not signifi- 
cantly affected by the presence of the Rydberg atoms. We 
present measurements on the excitation dynamics which 
show a significantly smaller Rydberg fraction in the con- 
densed sample than in the thermal sample. A simplified 
superatom model reproduces the main features of the 
measurement as a consequence of the changing density 
distribution upon condensation. The coherence of the 
condensate was not relevant for the understanding of the 
measurement of the overall excitation dynamics. Future 
studies will focus on making use of the coherence prop- 
erties of the condensate, e.g. in order to measure the 
spatial density-density correlation function of the Ryd- 
berg atoms by means of matter wave interferometry. 
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